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ABSTRACT 


The thermal evolutions of the Moon/ Mars, Venus and 
Mercury are calculated theoretically starting from cosmo- 
chemical condensation models. An assortment of geological, 
geochemical and geophysical data are used to constrain both 
the present day temperatures and the thermal histories of 
the planets' interiors. Such data imply that the planets 
were heated during or shortly after formation and that all 
the terrestrial planets started their differentiations early 
in their history. 

Initial temperatures and core formation play the most 
important roles in the early differentiation. The size 
of the planet is the primary factor in determining its present 
day thermal state. The Moon, smallest in size, is character- 
ized as a differentiated body with a crust, a thick solid 
mantle and an interior region which may be partially molten. 

It is presently cooling rapidly and is relatively inactive 
tectonically. 

Mercury, which probably has a large core, may have a 
500 km thick solid lithosphere and a partially molten core 
if it is assumed that some heat sources exist in the core. 

If this is not the case, the planet's interior temperatures 
are everywhere below the melting curve for iron. The 
thermal evolution is dominated by the core separation and the 
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high conductivity of iron which makes up the bulk of Mercury, 

Mars, intermediate in size, is assumed to have differ- 
entiated an Fe-FeS core. While the formation of an early 
crust is evident, large scale melting and differentiation 
of the mantle silicates has occurred at least up until 1 
billion years ago. Present day temperature profiles indicate 
moderate tectonic activity at the present time. 

Venus is characterized as a planet not unlike the Earth 
in many respects. Core formation has occurred probably during 
the first billion years after the formation. The present 
day temperatures indicate a partially molten upper mantle 
overlain by a 100 km thick lithosphere and a molten Fe-Ni 
core. We can expect that today Venus may have tectonic 
processes similar to the Earth's. 


1 . 


1 . INTRODUCTION 

The evolution of a planetary body is very strongly 
controlled by its thermal regime. Differentiation, volcanism, 
tectonic activity, and even magnetic field are all processes 
controlled by the internal temperatures and the thermal 
history. The calculation of the ' thermal evolution, however, 
requires strong constraints if it is to be realistic.. In 
recent years significant new data from lunar and planetary 
missions have become available. For the Moon these have 
provided constraints both on the present day temperatures 
and the past history of the thermal state of the lunar 
interior. For Mars, Venus and Mercury the thermal constraints 
are primarily qualitative. The purpose of this paper is 
to review these data and the thermal evolution models of 
the Moon and the terrestrial planets. 

An extensive number of studies on thermal history of 
the Moon and planets have been carried out (Urey, 1951, 1952, 
1955, 1957, 1962? Lubimova, 1958? MacDonald, 1959, 1963; Levin, 
1962? Anderson and Phinney, 1967; Fricker et al . , 1967, 1974; 
McConnell et ai. , 1967; Lee, 1968; Hanks and Anderson, 1969; 
Majeva, 1969? Urey and MacDonald, 1971; Papanastassiou and 
Wasserburg, 1971a; Hays, 1972; Wood, 1972; Reynolds et al. , 

1972? McConnell and Gast, 1972? Toksdz et al. , 1972c?- Tozer, 
1972a, 1972b, 1973? Hanks and Anderson, 1972? Duba and Ringwood, 
1973? Cassen and Reynolds, 1973? Toksdz and Solomon, 1973; 

Toksttz and Johnston, 1974? Johnston et al. , 1974; Siegfried 
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and Solomon, 1974; Schubert and Young, 1974), 

In the early studies the Moon and the planets were 
generally modeled with chondritic compositions and radio- 
activity. with the successful lunar missions and new 
geochemical data, it became necessary to revise the constraints 
used in lunar thermal calculations. The models generated 
since 1971 have incorporated these new data and generally 
give similar results. Generally, calculations for the other 
planets reflect this new knowledge about the chemistry of the 
solar system. 

In this paper then, we present thermal evolution models 
for the Moon, Mars, Venus and Mercury for a comparative 
study. We use similar computational techniques and thermal 
factors associated with melting, differentiation and core 
formation. These are discussed in Section 2. The thermal 
evolution of the Moon with the appropriate constraints are 
given in Section 3 and those of the planets in Section 4 of 
this paper. 

2. COMPUTATIONAL PROCEDURE AND MODEL PARAMETERS 

For the theoretical calculations of thermal evolution 
in a planetary body it is necessary to specify the initial 
conditions, heat sources, heat transport mechanisms, and a 
number of other parameters such as specific heat, density, 
and melting temperatures as well as the variation of these 
with temperature and pressure inside the planet. Generally, 
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not all of these parameters can be determined and theoretical 
calculations cannot define thermal models uniquely. 

The first input in theoretical calculations is the initial 
conditions assumed. Included in this is the heating due to 
accretion/ tidal forces/ induced electric currents/ and short- 
lived radioactivity, although the relative importance of each 
may vary. Next is the definition of the long term heat sources 
that are important throughout the age of the planet. These 
include heating due to the decay of radioactive isotopes and 
the possibility of long term heating due to tidal dissipation 

* 

and electrical currents. Finally, using the above inputs com- 
bined with other physical parameters, the conservation of 
energy equation must be solved numerically with the addition 
of certain perturbations designed to simulate melting, differ- 
entiation and convection. 

In this section we first describe the input parameters 
and then discuss the computational procedures involved in 
solving the thermal evolution problem. 

2.1 Model Parameters 

2.1.1 Initial Temperature 

The initial temperature of a planet is both a function of 
its mode of formation and of its immediate environment during 
and shortly after its origin. As will be discussed in Sec- 
tions 3 and 4, constraints on lunar and planetary thermal evo- 
lution favor initially hot and possibly molten models.. This 
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is especially true for the Moon. Energy sources which might 
have provided the heat to raise the temperature of the planets 
to solidus or near solidus levels include gravitational energy 
due to accretion, tidal dissipation, solar wind flux, short- 
lived radioactivity, and heating due to adiabatic compression. 
The temperature rise due to adiabatic compression amounts only 
to a few tens of degrees in the Moon but increases with the 
increasing size of the planet. The effect of short-lived 
radioactivity depends on the time between nucleogenesis and 

lunar formation. Also, Schramm et al. (1970) have shown that 
26 

Al , the most likely heat source in the early history of the 
solar system (Fish et al. , I960), had little effect on the 
initial heating of the Moon and meteorites. 

Tidal dissipation and solar wind flux could have been 
effective under special circumstances. The amount and the dis- 
tribution of heat is difficult to quantify in both cases. For 
the Moon, as an example, the effects of tidal dissipation would 
be small (Kaula, 1966), unless the Moon’s orbit was originally 
much closer to the Earth. This might have occurred shortly 
after a capture event or after the Moon formed in close prox- 
imity to the Earth (Singer, 1970). Hallam (1973) has calcu- 
lated that tidal dissipation, even with Q = 10, raises the 
internal temperature by only 50-100°K in 1.5 b.y. Even if 
the Moon was once close to the Roche limit, it would recede so 
rapidly that the effect of tidal heating would be negligible. 
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Heating by electric currents produced by a unipolar gene- 
rator driven by the solar wind might have been significant if 
the Sun passed through a T-Tauri stage (Sonnett and Colburn, 
1970 ) . • This effect is difficult to calculate and include 
explicitly in any thermal model. However, combinations of 
all the above effects may contribute to the initial temperatures 
We considered these effects by means of a "base" temperature 
added to accretional thermal profiles discussed below. 

The most promising initial heat source for the 
terrestrial planets is the gravitational energy of accretion. 
Classically this effect is described by the equilibrium 
of added gravitational potential energy with black body 
radiation and internal heat: 


GMfr) 44 

P — ■ £ ■■ ■ dr » ea(T(r) -T*) dt + pC p (T(r)-T b ) dr (1) 

where t is time, p is the density of the accreting particles, 

G is the gravitational constant, M(r) is the mass within 

radius r, a is the Stefan-Boltzmann constant, C is the 

P 

specific heat, and e is the emissivity (taken to equal 1) . 
T(r) is the temperature at radius r and T. is the base 

D 

temperature or the temperature of the accreting particles, 
including any effects of short-lived radioactivity. By 
specifying the accretion rate dr/dt and the physical 
parameters p, C , and , temperature versus depth for an 
accreted body c an be calculated from equation (1) . An 
accretion rate suggested by Hanks and Anderson (1969) and 
Mizutani et al. (1972) is used in this paper and is given by: 
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~ = ct 2 sin y t (2) 

dt 

The constants c and y are determined from the planet's final 
radius and the duration of accretion. Although the accretion 
rate is difficult to establish, it is a fair statement of the 
presumption that the rate grows as the accreted body grows and 
then tapers off as primary material is depleted. 

The accretion time is considered arbitrary, except that 
it is assumed to-be longer for larger planets. Generally the 
accretion time and T fa are adjusted to provide an initial tem- 
perature profile consistent with the constraints that require 
initial differentiation. For the Moon, it can be shown from 
equations (1) and (2) that the temperature will never reach 
the basalt solidus for accretion times greater than 1000 years, 
unless the emissivity of the lunar surface is decreased. This 
conclusion was reached independently by Mizutani et al. (1972) 
from a more detailed study of the accretion process. 

2.1.2 Heat Sources 

The most important long term heat source in thermal his- 
tory calculations is the heat generated from the decay of the 
long half-life radioactive isotopes U 238 , U 235 , Th 232 , and K 40 . 
Abundances of these isotopes have been measured in many Earth 
rocks and returned lunar samples (LSPET, 1972, 1973; O' Kelley 
et al. , 1970; Clark and Keith, 1973; Silver, 1973; Tatsumoto 
et al., 1973? Perkins et al., 1970; and others), and indirectly 
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by lunar orbital y-ray spectrometers (Metzger et al. , 1972). 
While surface concentrations vary greatly, we may conclude that 
these values must be much higher than those in the interior 
of the Earth or the Moon (Papanastassiou and Wasserburg, 1971b; 
and others) . The distribution of the isotopes in space and 
time are difficult to constrain, however, due to complex 
melting and crystal-liquid fractionation processes. 

Figure 1 shows total K and U abundances measured for 
Apollo samples (including crystalline rocks, soils, and 
breccias), some meteorites and terrestrial rocks. Figure 2 
shows Th/U versus U for the same suite of measurements. Al- 
though there are small and systematic deviations, the K/U 
ratio for lunar material averages around 2000, compared with 
10,000 for terrestrial rocks and 80,000 for chondrites. The 
low lunar K/U ratio is consistent with general depletion of 
volatile elements (LSPET, 1969, 1972) . From this, it is 
clear that neither terrestrial nor chondritic abundances can 
be used as the bulk radioactivity for the Moon. Furthermore, 
using Rb abundances for a model of the primitive lunar 
material (Papanastassiou and Wasserburg, 1971a) we can also 
reject chondrites. The Th/U ratio is consistently between 3.0 
and 4.0, for terrestrial, lunar and meteoritic material. 

It is essential then, to determine or estimate the bulk 
radioactive concentrations for each of the planets. In the 
case of the Moon, given that chondrites and terrestrial rocks 
are not good approximations of the bulk radioactivity, we 
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turn to achondrites - specifically howardites and eucrites. A 
q ene tic relationship between achondrites and the Moon has 
been suggested by several investigators (Duke and Silver, 

1967; Anderson, 1973; W&nke et, al_., 1973). Howardites have, 
on the average, a uranium concentration of 35 ppb with K/U 
= 9800, while brecciated eucrites, enriched in refractory 
elements have U = 130 ppb with K/U - 2300 (Mason, 1971) . 

It is interesting to note that the K/U ratio for eucrites 
is close to that of the lunar material while the ratio fur 
howardites is somewhat higher. 

If we take K/U for the bulk radioactive concentration 

of the Moon to be 2000, and Th/U to be 3.5, then the average 

radioactivity of the Moon is a function of the bulk uranium 

content. Given the measured Apollo 15 and 17 heat flows of 

2 

about 30 ergs/cm -sec (Langseth et al. , 1972, 1973) , Toks&z 
and Solomon (1973) have shown that the bulk uranium concen- 
tration must be 60 + 15 ppb to match the heat flow value. 

In this paper, then, we let U = 60 ppb with K/U » 2000 and 
Th/U = 3,5 for the Moon. 

For Mars, Venus and Mercury we will utilize values based 
on condensation models of the solar nebula (Lewis, 1972; 
Grossman and Larimer, 1974; Anderson, 1973). These are shown 
in Table 1. In the calculations we consider variations for 
individual cases as discussed in later sections. However, 
expect that for Venus, terrestrial abundances are a good 


s 


we 
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approximation. For Mars we use values between terrestrial and 
chondrite abundances (Johnston et al . , 1974).’ For Mercury, 
the uranium values are close to those of the Moon but the 
planet is most likely lacking in K 40 (Siegfried and Solomon, 
1974). The relevant decay constants, decay energies and 
isotopic abundances are taken from Clark (1966) 

2 -l*3 Distribution of Heat Sources 
In most cases, the thermal models presented in this 
paper assume an initially homogeneous distribution of heat 
sources. However, the constraints on the internal temperatures 
rule out at least for the Moon a present day uniform distri- 
bution of heat sources. When melting occurs, one would 
expect the magma to be enriched in U, Th, and K and to 
subsequently transfer these heat sources to the surface. To 
account for this, at discrete time steps the heat sources from 
the molten zones are differentiated upwards, leaving primordial 
radioactive abundances in any solid region below the zone. 

Heat sources differentiated towards the surface are allowed 
to decrease exponentially with depth so that: 

H(r) = A q exp -iSgElj (3) 

where A Q is the surface abundance, R is the radius and H(r) is 
the radioactive abundance at radius r. Such a distribution 
has been formulated for the Earth (Lachenbruch, 1968). The skin 
•depth h decreases with time and is adjusted so that the present 
day surface radioactivity, if known, is approximated. 
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2.1.4 Thermal Conductivity and Other Parameters 
In a differentiating planet, heat is transported by con- 
duction and by convection. At temperatures far below solidus, 
the role of convective flow is small, and heat transport is 
governed by thermal conductivity. This quantity is tempera- 
ture dependent and is taken to be the sum of two terms, lattice 
(or phonon) conductivity and radiative conductivity. The 
relative dominance of each term is of crucial importance in 
solving the conduction problem of the planets. Although effect- 
ive thermal conductivities have been measured for terrestrial 
and lunar materials (Figure 3), these do not always represent 
the bulk properties of the Moon of planetary interiors. A 
theoretical formulation of conductivity as a function of tempera- 
ture by MacDonald (1963) assumed the opacity to be constant, 
resulting in a cubic temperature dependence of the radiative 
term and thus its dominance over all temperatures. However, 
Schatz and Simmons (1972) proposed that for a hypothetical dun- 
ite, radiative transfer is suppressed by grain-boundary ex- 
tinction, resulting in an increase in the opacity with tempera- 
ture. At low temperatures, then, tne effective conductivity is 
dominated by the lattice term which decreases with increasing 
temperature. At higher temperatures the radiative term, in- 
creasing roughly linearly with temperature, becomes dominant. 

This model has been used in previous thermal history calculations 
at M.I.t. (Toksdz et al . , 1972c;Toksdz and Solomon, 1973). 
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However, for the Moon, upon examination of Figure 3 and values 
of conductivity in many previous lunar models, one could just 
as well take the effective conductivity to be a constant. 

In some of the lunar models and for Mercury's mantle we shall 
treat it thus, with K = .45 x 10 6 ergs/cm-sec °C while other 
calculations assume the model of Schatz and Simmons (1972). 

The total conductivity K is given by 


K = K_ + K_ 

L R 

where — lattice conductivity and ~ radiative conductivity 
As functions of temperature T, these are given by 

K = 4 - 184 x 10 7 
L 30.6 + .2 IT 


and 


K 


R 


0 for T f 500 °K 
230(T-500) for T > 500°K 


Other physical parameters required for the calculations 
are the melting curve, surface temperature, specific heat, heat 
of fusion, and density. A summary of the parameters common to 
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all models is given in Table 1. Appropriate melting curves 
are used for each planet taking into account the composition 
and pressure. 

2.2 Computational Techniques 

In this section we describe the thermal calculations for 
a spherically symmetric planet, taking into account conduction, 
melting, simulated convection of molten material and differen- 
tiation. 

2.2.1 Conductive Models 

Assuming that conduction is the only mode of heat trans- 
fer in solid material (an approximation to convection 
will be discussed later) , thermal evolution models are 
calculated using the finite difference solution of the heat 
conduction equation: 

e c P H " 4 IF { * 2 « ^ + «<*'*> < 4 > 

r r 

where is the specific heat, T is the temperature, r is 

the radius, K is the thermal conductivity and H(r,t) is the 

♦ 

heat source term. The latter two quantities were discussed 
above . 

The details of the calculations are described fully in 
the appendix of ToksdSz et al. (1972c) and are stated here 
briefly. The finite difference analog of equation (4) which 
conserves heat flux is 
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where p - Ar, r t- np f t - mAt, a = At/pC and T™ = T(np, mAt) . 

r 

The stability condition relating the time increment and the 
grid spacing is: 

2 2 
n 

( 6 ) 


G P Ar 

At < 


2 K , ,L2 

max (n+~) 


where is the maximum value of the conductivity at a 

given time step. At is computed at each iteration using 
equation (6) with a factor of 4 instead of 2 in the denomin- 
ator. A grid spacing, Ar, of 20 km is used. 

Boundary conditions for the problem require the surface 
temperature of the planet to be constant, and = T™ where 
T 1 is the temperature at r x =» Ar and T™ is the temperature 
at r 2 = r l + * The solut ion starts at r = Ar to avoid 

a singularity at r = 0. Surface heat flux is calculated 
by adding the heat flow at the boundary between the upper 
two finite difference cells to the steady state heat flow 
due to heat production in the uppermost half cell. 

2.2.2 Simulation of Melting and Convection 

The effects of melting and convection of molten material 
on temperature within the body can be modelled following the 
technique described by Reynolds et al. (1966) . In this 
method the temperature at each grid point at each time step 
is compared to the melting temperature at the respective 
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grid point. If the temperature is greater than the melting 
point, the difference divided by the density and specific 
heat gives the heat equivalent, AH. If AH is less than or 
equal to the heat of fusion, AH^, the material at that grid 
point is considered partially molten, and the temperature 
is set equal to the melting point. This process is continued at 
each iteration until the accumulated AH is greater than AH^. 

At this point in time the material is considered to be 
completely molten and the temperature is allowed to rise 
above the solidus. The same process allows completely 
molten material to solidify releasing heat equal to AH^. 

In simulating convection of molten material, the same 

technique as outlined above is followed, with the exception 

that the temperature of the completely molten material is 

constrained to the melting curve. Any increase in temperature 

that would raise the temperature above the melting point 

is converted to its heat equivalent and transferred upward 

to the next grid point and is considered to be a temperature 

rise due to convection. The conversion to heat equivalents 

is necessary due to the increase in the volume element with 

increasing radius. The advantage of this approximation to 

* 

convection is that it does not require the inclusion of mass 
transfer terms in the heat conduction equation. This scheme 
is not limited to convection of completely molten material. 
Convection of partially molten material is easily approxi- 
mated by reducing the heat of fusion. Convection by solid 
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state creep requires techniques different from these described 
above. Various attempts at approximating solid state convec- 
tion and determining its relative effect on the thermal evolu- 
tion of the Moon and planets are discussed in the following 
section. 

2 . 3 Solid State Convection 

The question of whether or not thermal models calculated 
assuming conduction mechanisms are unstable to solid state 
convective motions has been discussed in many papers, most 
recently by Tozer (1972a, 1972b), Turcotte et al. (1972), 

Cassen and Reynolds (1973a, 1973b), and Schubert et al. (1974). 
Much of this work has been notivated by Schubert et al. (1969) . 
Using a number of conduction temperature profiles for the 
planets, these investigators applied a linear stability analy- 
sis and concluded that thermal convection would occur in each 
case. Earlier, Runcorn (1962, 1967) attempted to explain de- 
partures of the lunar surface from hydrostatic equilibrium 
by employing convection. Whether or not convection has played 
a dominant role in the evolution of the planets, however, is 
still debatable. The question is not easily answered since 
the presence (or absence) of convection is a complex function 
of temperature and the rheology of the planetary material. 
Furthermore, in the case of the Moon, the characteristics of 
the mare basins and their possible filling with basalts by gra- 
vitational upwelling of magma, implies the presence of more 
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complicated geometries than the simple cellular convection 
pattern postulated by some investigators. Such 
upwelling type heat transport is well approximated in our 
calculations, described in Section 2.2.2. Also, the attenua- 
tion of shear waves and the presence of deep moonquakes in the 
lunar interior as discussed in Section 3.1.3 present diffi- 
culties for those who hypothesize solid state convection in 
the Moon. 

The strongest advocate for the dominance of solid state 
creep in the Moon and planets has been Tozer {1972a, 1972b, 
1973). Tozer has, in all his papers, postulated that convec- 
tion is at least as important as conduction for any planetary 
body of radius greater than a few hundred kilometers. This 
theory, based on the assumption that at any finite temperature 
an arbitrarily small non-hydrostatic stress produces some per- 
manent deformation of the material, suggests that at some vis- 
cosity, or range of viscosities, this motion, or rate of creep, 
becomes "fast" enough to provide a very efficient method of 
heat transfer. 

In this paper we do not go into detailed convection calcu- 
lations. In some cases, however, we include models that approx- 
imate solid state creep by arbitrarily raising the thermal con- 
ductivity of material with temperatures above a certain thresh- 
old value. 
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3. THERMAL EVOLUTION OF THE MOON 

3 . 1 Constraints on Thermal Evolution 

Thermal evolution calculations generally require boun- 
dary conditions which would indicate directly or indirectly 

the temperatures at certain times inside the planetary body. 
In addition to these, data are needed to specify material 
behavior, heat transport mechanisms and heat sources. 
Normally, insufficient data and constraints are available. 

Thus calculations generally are carried out with a set of 
initial conditions and assumptions, and the conditions are 
modified if the temperature models fail to satisfy available 
constraints. 

Data and samples from Apollo missions imposed some 
strong constraints on the thermal evolution of the Moon. Some 
of these such as heat flow, electrical conductivity, seismic 
velocities, attenuation and tectonism, pertain to present- 
day temperatures. Some others such as the presence of a 
^ enti ated lunar crust and chronology of lunar igneous 
and magmatic activity indicate the thermal state early in the 
lunar history. 

3.1.1 Lunar Heat Flow 

Two heat flow measurements have been made on the Moon, 
one each at the Apollo 15 (Hadley Rille) and Apollo 17 
(Taurus Littrow) sites. The experiment and the results are 
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described by Langseth et al_. (1973). The heat fluxes have 
been determined by measuring temperature gradients and ther- 
mal conductivities in bore holes to a maximum depth of 
234 cm. In that depth range temperature gradients in the 

lunar regolith are between 1.3 and 1.7 °K/m, and the thermal 

, -4 -4 

conductivity is m the range 1.7 x 10 and 2.0 x 10 W/cm-°K. 

The heat flow through the surface based on the measured 

2 

temperature gradients and conductivity is 31 ergs/cm sec 

2 

at the Apollo 15 site and 28 ergs/cm sec at the Apollo 17 
site (Langseth et al . , 1973) . In each case the estimated 
uncertainty is about ±20%. 

The Taurus Littrow and Hadley Rille sites are located 

in the basalt-flooded valleys at the rims of the Serenitatis 

and Imbrium basins, respectively. There are no direct 

data to suggest that the thermal characteristics of these 

sites are any different from other parts of the Moon. Thus 

an average heat flow value of about 30 ergs/cm 2 sec will be 

adopted for the Moon in these calculations . This value 

is about one-half of the average heat flow of the Earth 
2 

(63 ergs/cm sec) . If we assume steady state conditions for 
both bodies and take into account the differences in the radii 
(average radii: ^ oon = 1738 km, R Earth = 6371 km), surface 

area (S Moon /S Earth = 0.074) and the volume ■(y |loon /V Earth = 0.02) 
we find that the lunar heat flow is higher by about a factor 
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of 2 compared to the Earth. As was discussed earlier in 
the article, this implies higher heat generation and radio- 
active sources inside the Moon than inside the Earth. 

Another important measurement arising from the heat 
flow experiments is the value of the near surface temperature. 
At a depth of about 70 cm, the measured average ambient 
lunar Temperature is 253°K. This will be used as the surface 
temperature in our calculations. 

3»1»2 Electrical. Conductivity and inferred Temperatures 

Three magnetometers deployed on the lunar surface, one 
each at the Apollo 12, 15 and 16 sites, have recorded the 
magnetic fields induced in the Moon by large scale transient 
events. In addition to these, the ambient and time-dependent 
magnetic fields in the lunar environment are measured by 
the Explorer 35 satellite magnetometer which is orbiting the 
Moon. 

The analysis of the transient fields recorded simul- 
taneously by the Explorer 35 magnetometer and the surface 
magnetometer has made it possible to calculate the electrical 
conductivity inside the Moon (Dyal et al . , 1972, 1973; 

Sonett et al . , 1971, 1972). The Moon responds to transient 
magnetic events which induce eddy currents in the lunar 
interior. The response or transfer function is 
computed from the induced field at the Moon's surface and 
the forcing (source) function recorded by the Explorer 35 


magnetometer . 
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These results are * interpreted by calculating the 
theoretical response of the Moon with radially varying 
conductivity. Although there are some limitations of the 
theoretical calculations and some difficulty in matching 
all the available data, the results definitely show an 
■increase of electrical conductivity with depth into the Moon. 

To a depth of about 1000 km the conductivity profile can be 
characterized by three layers (Dyal et al. , 1973) or by a 
greater variety of models including current layers (Sonett 
et al. / 197 2) . The resolving power of the data is not 
sufficient to detail the conductivity below this depth. 

Conductivity depth profiles given by Dyal et al. (1973) 
and the three layer model (a good representation for the 
average models) of Sonett et al. (1972) are shown in 
Figure 4. Basically the models require very low conductivity 
- 10 mho/m) for the outer 50 km of the Moon. For the 
intermediate layer (1100 < R < 1700 km; R is radius) the 
conductivity is in the range of a 2 = 1 to 7 x 10~ 4 mho/m. 

For the deeper interior there is a greater discrepancy* 

While Sonett et. al. (1972) give a value of about - 10~ 3 mho/m 
for R < 1200 km from the sun-lit side data, Dyal et a_l. (1973) 
give o 3 - 10~ 2 for R < 1000 km. 

From the conductivity profiles the temperatures inside 
the Moon can be estimated if the composition and temperature 
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dependence of conductivity are known. For a given composition, 

, A . 

conductivity, increases with temperature: - o Q 

where A = activation energy, k = Boltzmann constant, T = 
absolute temperature , o Q is a constant dependent on the 
material. If olivine were assumed for a composition, and 
a lower bound on the conductivity of olivine is used, a 
temperature of about 1000°C is estimated for the Moon, at 
about 700 km depth (Dyal et al . , 1973 ; Sonett et al_. , 1972). 

If newer laboratory data on pyroxenes and olivines - the 
most likely candidates for mantle composition (Duba and 
Ringwood, 1973; Housley et al . , 1972; Olhoeft et al . , 1973) 
are used, higher temperatures (i.e. 1000°C at 400 km and 
1400° at 800 km depth) are estimated (Anderson and Hanks, 1972; 
Duba and Ringwood, 197 3) . Some bounds on temperature models 
are shown in Figure 4; in all cases temperatures are below 
the solidus to a depth of about 600-700 km. At greater 
depths, temperatures either approach or exceed the solidus, 
or some changes in the composition increase the conductivity. 
Another feature of this deep region is that both conductivity 
and temperature profiles flatten. 

3.1.3 Structure of the Lunar Interior and Seismic 
Wave Attenuation 

Seismic studies of the Moon have revealed much 
valuable information on the structure and the rheological 




22 


properties of the lunar interior. The presence of a 
layered structure indicates differentiation, and melting in 
the past. The attenuation characteristics of seismic waves 
may be explained in terms of increased temperatures or 
partial melting at the present. 

The seismic velocity structure for the outer 150 km of 
the Moon has been determined from the analysis of man-made 
impact seismograms recorded by the Lunar Seismic Network 
(Toksbz et al . , 1972a, b, 1973a). This structure shows the 
presence of a lunar crust overlying a mantle. Over the mare 
the crust may be two-layered (the upper layer being mare 
basalt) with a total thickness of about 60 km. In the 
highlands and on the backside, the crust may be somewhat 
thicker. The compositional inferences that can be made 
from the velocities and other data suggest that the crust 
is chemically different from the mantle and it may have 
formed by the upward differentiation of the lighter minerals. 

Below the crust the lunar mantle may extend to a 
depth of about 1000 km, probably with nearly constant com- 
pressional velocities (Nakamura et_ al. , 1974; Toksdz et^ al. , 
1974) . Data from moonquakes and meteorite impacts specify 
the properties of the mantle. The recording of well-defined 
shear waves from the deep moonquakes implies that the mantle 
must be sufficiently rigid to 700 or 1000 km depth to prevent 
appreciable shear wave attenuation. 
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Seismic shear waves from farside impacts and farside 
moonquakes that have penetrated deeper than about 700 km 
into the Moon are attenuated (Latham et al . , 1973; Nakamura 
et al . , 1973; Toksftz et al . , 1974) . This indicates a 
possible "softening" of the material in the central region 
of the Moon below about 700 or 1000 km. Any softening which 
would reduce Q to less than about 500 would explain the 
observations. This could be achieved by temperatures approaching 
the solidus, by a very small amount of partial melt or by 
other mechanisms (such as perhaps an increase in the amount 
of volatiles in the deep lunar interior or a different bulk 
composition) . 

If we assume that the attenuation is due to temperature 
effects, this would place some limits, within the compositional 
constraints, on the temperatures below 700 km in the Moon. 

For a pyroxene-olivine rich composition this is 1600 °C. If 
the Moon has an Fe-FeS rich core (see Solomon, 1973 for a 
discussion) the minimum temperature for it would be about 
1000 °C given by the Fe— FeS eutectic. Thus the present-day 
lunar "core" temperatures are between about 1000-16 00 °C. 

3*1.4 Constraints Implied by Lunar Viscosity, Seismicity 
and Mascons 

The lunar interior described in the previous section 
and shown in Figure 5 can be divided into three units: a 
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crust, a relatively rigid 1000 km thick "lithosphere" and a 
relatively soft deep interior "asthenosphere" . 

The crust and the lithosphere have supported lunar 

mass anomalies ( "mascons" ) for over three billion years. 

This requires that (1) at the time of the raascon formation 

and mare flooding, the lunar "lithosphere" was thick enough 

to support the loads (Solomon and Toksttz , 1973) and (2) 

the viscosity of the "lithosphere" must be greater than about 
2 6 

10 poise (Baldwin, 1970; Arkani-Hamed , 1973). This high 
viscosity requires temperatures significantly lower than the 
melting curve in the outer few hundred kilometers of the 
Moon . 

Another important constraint on the present-day thermal 
state is that the lunar interior is the very low level of 
seismic activity in the Moon. Total seismic energy release 
in the Moon is about 10 11 ergs/year, thirteen orders of 
magnitude lower than that of the Earth (Latham et ad., 1972, 
1973; Lammlein et ad., 1974). All moOnquakes are very small 
and those whose focal depths can be determined have hypo- 
centers between 700 and 1200 km depth. These indicate 
that if there is any convection or tectonic motion taking 
place in the Moon today, this must be below about 700-1000 km. 
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3.1.5 Chronology of Lunar Igneous Activity 

The time history of differentiation and igneous activity 
imposes the strongest constraints on the early thermal state of 
the Moon. The chronology of the lunar igneous activity starts 
f rom ■ the formation of the original crust about 4.6 billion 
years (b.y.) ago as deduced from model ages of lunar rocks and 
soils (Silver,. 1970; Wasserburg and Papanastassiou , 1971; Tat- 
sumoto £t ?il . , 1972; Tera et al . , 1973). More recently a Rb-Sr 
age of 4.6 b.y. has been obtained for an Apollo 17 dunite 
fragment (72415) by Albee et al . (1974) . To create a 60 km 
thick feldspathic crust described in the previous section re- 
quires the total or partial melting and differentiation of at 
least one-half of the total volume of the Moon (Wood et al. , 
1970). This requirement strongly constrains the initial tem- 
perature of the lunar interior. 

Following the formation of the initial crust, the 
chronology of lunar rocks (an excellent summary has been 
compiled by Doe, 1973) is defined by three episodes in the 
lunar evolution: 

1- The period 4.4 - 4.0 b.y. ago was characterized by 
continuing igneous activity, brecciation by meteoritic 
bombardments and metamorphism of the older rocks. Because 
of the subsequent events, the detailed record of this 
episode has been largely obscured. Available data include 
ages of anorthosites (4.1 b.y. for lunar rock 15415, Husain 
et al., 1972a), the formation ages of KREEP basalts 
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(Papanastassiou and Wasserburg, 1971b; Schonfeld and 
Meyer, 1972) and. a number of the "recrystallization" ages 
of the metamorphosed breccias from the Apollo 12, 14, 16, 

17 and Luna 20 sites (See Doe, 1973 for a compilation). 

These are shown in Figure 6. 

2. Episode of basin excavation (4.0 - 3.9 b.y. ago). 

This period is characterized by intense bombardment, re- 
setting of the geologic clocks, and the excavation of the 
major basins such as Imbrium, Tranquilitatis , Serenitatis, 
Nectaris, etc. ( Tera et al . , 1973? Tera and Wasserburg, 1973). 

3. Episode of mare flooding (3.9 - 3.16 b.y. ago). 

The basins formed were filled with mare basalts during this 
period. The ages of these lunar basalts are well documented 
by dating samples from the Apollo 11, 12, 15, 17 and Luna 

16 missions (Wasserburg and Papanastassiou, 1971; 
Papanastassiou, 1971a, 1972? Tatsumoto et al. , 1972? 

Murthy et ad., 1973; Turner et al_. , 1973; Husain et al. , 

1972b, and others) . 

The failure to discover lunar igneous rocks younger 
than about 3.16 b.y. is an extremely important constraint 
in the lunar thermal history. Since basalts were probably 
derived from sources to about 400 km depth, the termination 
of lunar magmatic activity represents the end of melting 
or partial melting in the upper mantle of the Moon. 



During the past 3.0 b.y. , the lunar surface has 
continued to evolve as a result of meteorite impacts and 
cratering . There have been local melting and metamorphism 
in this period. There is no direct evidence for extensive 
volcanic activity or basalt flooding in the past 30. b.y. 

3.1.6 Magnetism of Lunar Rocks 

Although the Moon does not have a measurable dipole 
magnetic field at the present, the lunar crust is extensively 
magnetized. This remanent magnetization has been verified 
by magnetic field measurements from the orbiting Apollo 
sub- satellites (Coleman et al . , 1972) , on the surface 
(Dyal et al., 1972), and from the remanent magnetization of 
returned lunar samples (See Fuller, 1974; Nagata et al. , 1972; 
Pearce et al . , 1972; Collinson et al . , 1972, and others). 

This stable remanent’ magnetization was most likely acquired 
on the Moon as the rocks cooled through the Curie point of 
iron in the presence of a magnetic field. The intensity of 
the ambient magnetic field at the time of the formation of 
these rocks is estimated to be greater than 1000 gammas, 
and probably as high as a few thousand gammas (Gose et a_l. , 
1973) . 

The origin of the ancient magnetic fields, which 
magnetized the initial lunar crust and the lunar rocks, is one 
of the most important and as yet unresolved problems in lunar 
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studies. From the ages of the crust and the rocks, the 
magnetizing field must have been present at least from 4.6 
to 3.2 b.y. ago. The hypothesis of an external 
magnetizing field is not favored because of the difficulties 
associated with having a steady field near the Moon of 
1000 gammas or more for a period of 1.5 b.y. 

The favored hypothesis of the internally generated 
magnetizing field would require either a lunar dynamo or 
magnetized lunar interior (Runcorn and Urey, 1973; Runcorn 
et al . , 1971; Strangway et. ad., 1973) . In the dynamo model, 
it is necessary for the Moon to have a conducting iron or 
Fe/FeS core of sufficient size early in its history. The 
formation of such a molten core requires high internal 
temperatures early in the lunar history. In the case of a 
highly magnetized model of the lunar interior (Runcorn and 
Urey, 1973; Strangway et al . , 1973) it is required that the 
temperatures below a few hundred km depth remain below the 
Curie temperature (about 760°C) while all the magmatic 
activity and differentiation take place at shallower depths. 
Such a model would put severe requirements on the thermal 
evolution models and also require a high magnetizing field 
at the time of the lunar accretion. 

Thus the magnetic history of the Moon according to the 
above models would require either extensive melting and dif- 
ferentiation of the core at the onset of lunar evolution, or 
temperatures less than about 800°C for the first 1.5 b.y. in 
the lunar mantle. 
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3.2 Temperature Models for the Moon 

In this section we present several specific models 
for the thermal evolution of the Moon which illustrate the 
effect of varying initial conditions and bulk radioactivities. 
Previous papers (Toksftz e_t ad . , 1972c; Solomon and Toksttz, 
1973; Toksdz and Solomon, 1973) have discussed many of the 
models presented here and have assessed their ability to 
satisfy the constraints listed in Section 3.1. In this 
section we divide the models into three subgroups? initially 
homogeneous, convective and inhomogeneous models. 

3.2.1 Initially Homogeneous Models 

We begin by considering a rather unrealistic model: 
initially cold (T - 0°C everywhere) , that is a valuable 
exercise in the respect that it provides a verification of 
the computational scheme. Figure 7 shows three present day 
temperature profiles each for the given present day bulk 
uranium concentrations. Clearly, none of the models shown 
here satisfies the constraint of near surface melting in the 
first billion years of the Moon's history. In fact, the 
basalt solidus {dashed line) is not reached at any point in 
time in the lunar evolution. Furthermore, while reasonable 
uranium concentrations of greater than 37 ppb provide present 
day melting in the lunar interior, these models cannot 
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generate enough heat to satisfy the above constraint. Since 
radioactive decay alone cannot provide the heat necessary 
for early melting and differentiation, we must consider 
models with high initial temperatures. 

In Section 2 it was stated that one might consider the 

effects of early heating by tidal dissipation and/or solar 

wind flux as the Sun passed through a T-Tauri, stage by 

investigating an initially completely molten Moon. Such a 

model is shown in Figure 8. In this case, all the radioactive 

heat sources are immediately differentiated and concentrated 

towards the surface. The formation of a dense core (Solomon 

and Toksbz , 1973) composed of iron or FeS would easily be 

accomplished. Because of the depletion of heat sources in 

the lunar interior and the very effective heat transport in 

the convecting molten regions, this model cools rapidly and 

is completely solid after 2 b.y. At the present time it is 

completely solid and continuing to cool. With an average 

uranium composition of 60 ppb , the heat flow value is 
2 

29 erg/cm -sec, within the reported range of the Apollo 15 
and 17 measurements. Although this model satisfies the 
constraints of early initial melting and heat flow, it does 
not allow for a small partially molten core as implied by the 
seismic data. 

Since the normal process of accretion can produce initial 
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near surface partial melting without employing some 
catastrophic event, we now consider a model of lunar 
thermal evolution with initial temperatures calculated 
according to equations (1) and (2) . This then allows for 
early segregation of the lunar crust and differentiation of 
the radioactive heat sources. 

A model that satisfies all the major constraints presently 
known is shown in Figure 9. The initial temperature 
is calculated from the accretion model with an accretion time 
of 100 years and a base temperature of 80G°C, which allows 
for early heating by sources other than accretion. Initial 
melting extends to a depth of about 800 km, giving 
sufficient volume to differentiate a primordial lunar crust. 

The present day uranium concentration averages to 60 ppb . 

6 

The thermal conductivity is taken to be constant at .45 x 10 
erg/cm~sec-°C. 

It can be seen from Figure 9 that this model cools at 
the surface giving a solid lithosphere growing at a rate of 
about 220 kra per b.y. During the period of mare filling, 
the upper extent of melting progresses from 140 to 280 km, 
in agreement with the depth of origin of mare basalts (Ringwood 
and Essene, 1970) and with the necessity of maintaining a 
thick, cool lithosphere in order to sustain the stresses 
associated with mascori gravity anomalies. Shortly after 
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1 b.y. , a convecting core develops due to undifferentiated 
heat sources of nearly primordial radioactive abundances near 
the center of the Moon. The core has a maximum radius of 
1340 km at 1.5 b.y. and then slowly cools by conduction 
through the lithosphere. At present, this model is molten 
below 1000 km. This is in agreement with S v?ave attenuation 
in the core region and with the occurrence of moonquakes 

2 

at 700-800 km depth. The surface heat flow is 30 erg/cm -sec, 
consistent with the measured values. 

3.2.2 Convective Models 

We now consider an approximation to solid state 
convection in the Moon. We propose to simulate the efficient 
heat transport of convection by raising the thermal conduc- 
tivity by a factor of ten of a grid point that has a temper- 
ature greater than a critical temperature of 1000 °C. This 
is equivalent to Tozer’s steady state "core" temperature. 
Otherwise, the conductivity follows the scheme of Schatz 
and Simmons (1972). We permit this enhanced conduction after 
1 b.y., i.e. after rapid differentiation is completed and 
a lithosphere has formed. As in the case of convection of 
molten material, this approximation does not include mass 
transport phenomena, nor is it a statement on any particular 
solid state creep mechanism* It only enhances the heat 
transport. 
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Figure 10 illustrates a model with convection as 
described above. The initial conditions and present day 
bulk uranium concentration are the same as for the pre- 
vious model (Figure 9). The evolution proceeds as before 
until 1 b.y. At this point, convective processes dominate 
and the conductivity is raised for temperatures above 
1000°C. At the present, a convecting core of radius 1100 
km and temperature nearly 1000 °C underlies a nonconvecting, 
cool, rigid lithosphere. The present heat flow is 29 
erg/cra -sec. The features of this model are essentially the 
same as Tozer's models. 

3.2.3 Inhomogeneous and "Blanketed"' Models 

We have considered up until this point models which 
assume an initially homogeneous distribution of heat sources 
in the Moon. The only exception was that of the molten Moon 
in which all the radioactive isotopes were initially 
differentiated to the surface. One type of model which allows 
relatively cool interior temperatures and satisfies the 
heat flow constraint without invoking solid state convection 
has been suggested by several investigators (Gast, 1972; 

Green et al^. , 1972; and others) . This involves a decrease 
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in the concentration of radioactivity in the Moon with depth 
as a primary feature associated with inhomogeneous accretion. 

Figure 11 shows a thermal history model from Toksdz and 
Solomon (1973) similar in many respects to the one proposed 
by McConnell and Cast (1972) . The initial temperature 
profile assumes accretion in 1000 years with a base temper- 
ature of 0 °C , We assume that the radioactivity is stratified 
such that the present day uranium abundance is 120 ppb in 
the uppermost 200 km, 50 ppb in the next 250 km, and 10 ppb 
in the remainder of the Moon. This gives a bulk uranium 
concentration of 53 ppb. The K/U and Th/U ratios are fixed 
with values given in Table 1. The conductivity follows the 
Schatz and Simmons (1972) model. Although the initial 
temperature barely exceeds the solidus between 40 and 140 
km depth, shorter accretion times would provide initial 
melting at greater depths as in the previous models presented 
in this paper. With differentiation allowed for in four 
discrete time steps between 0 and .9 b.y., the evolution of 
the upper few hundred kilometers is similar to that in 

Figures 9 and 10 The present day heat flow is 32 erg/ 

2 

cm -sec and the interior of the Moon has temperatures well 
below the solidus. 

The difference between this model and the others 
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discussed earlier, is that the interior of the Moon is 
depleted in radioactive elements at cold temperatures, and 
thus heats slowly so that the Moon is everywhere solid at 
present. Whether or not a cold accretion can accompany 
primary chemical layering is debatable, however, such an 
origin is incompatible with a molten iron-rich core early in 
the Moon’s history. 'Furthermore, as in the case of solid 
state convection, the present day temperature profile does 
not seem to satisfy the constraint implied by the attenuation 
of S waves in the deep interior. 

We finally examine the effect of a low thermal conduc- 
tivity ’blanket” layer on the Moon. Measurements of the 
effective thermal conductivity of uncompressed lunar 
soils (Cremers, 1971) indicate values more than two orders 
of magnitude less than that assumed for the interior (see 
Figure 3) . Some investigators have argued that the widespread 
distribution of these lunar fines following the large 
impact events 4 b.y. ago might have produced localized 
melting at shallow depths. Figure 12 compares such a model 
with that of Figure 9, 100 m.y. after the soils are assumed 

to have been scattered. Since we are interested in the 
maximum possible effect of such an event, we will assume the 

lunar regolith to extend to a depth of 10 km with a thermal 

4 

conductivity of .45 x 10 erg/cm-sec-°C . Both values are 
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probably overestimates . The depth of dense fracturing in 
the lunar crust has been estimated to be less than 10 km by 
Toksttz et slI. (1973) . The thermal conductivity increases 
rapidly with the compression of the soils and in situ 
conductivity measurements made at the Apollo 15 and 17 heat 
flow sites (Langseth et al . , 1973) yield values only about 
one order of magnitude lower than the conductivity assumed 
for the bulk of the Moon. In any case, the maximum effect 
shown in Figure 12 amounts to a rise in temperature of only 
about 180 °C at 25 km. Thus it can be confidently concluded 
that a uniform thermal "blanket" would have had a minimal 
effect on the early evolution of the Moon. 

While the models discussed above are limited in number, 
they do serve to demonstrate the importance of initial 
conditions, geophysical constraints and heat transport 
mechanisms in determining the thermal history of the Moon. 

The range of present day temperature profiles inside 
the Moon which satisfy the main constraints listed in Section 
3.1 is shown in Figure 13* One encouraging aspect of this 
figure is that regardless of how the temperatures are calcu- 
lated, whether conduction or convection mechanisms are prevalent, 
the present day temperatures are relatively narrow. Hope- 
fully, with additional data, not only the present day 
temperature limits will be narrowed but also the early history 
of the Moon will be better understood. 
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3 . 3 Evolution of the Moon 

The thermal state of the lunar interior as a function 
of time is shown in Figure 14 along with major episodes in 
lunar evolution. During the first two billion years, the 
lunar upper mantle undergoes sufficient melting to account 
for the differentiation of the crust and the subsequent lunar 
vo lean ism and mare filling. The zone of melting deepens 
with time. The lithosphere thickens at the rate of about 
220 km/b.y. during this period. In particular, the 
shallowest melting progresses from 140 to 280 km depth during 
the period of mare filling, in agreement with the depth of 
origin of mare basalts (Ringwood and Essene, 1970) and 
with the need for a reasonably thick lithosphere to sustain 
the stresses associated with inascon gravity anomalies. The 
disappearance of melting in the mantle coincides roughly 
with the termination of magmatic events. At the present 
the whole Moon is cooling. The deep interior, below a depth 
of 1000 km may be hot enough for partial melting. This will 
explain the S-wave attenuation in this zone as described in 
Section 3.1. 

If indeed partially molten, the rheologic properties 
of the deep lunar interior may be similar to the Earth's 
asthenosphere and it may have convective motions. Such 
convection may exert small stresses at the bottom of the 
lithosphere where the moonquakes occur, but it could not 
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induce sufficient stresses in the thick, lunar 
"lithosphere'’ to cause large moonquakes or active tectonic 
motions. 

Present day temperatures shown in Figures 13 and 14 
generally exceed the melting temperatures of iron or iron/ 
iron sulfide combinations in the deep lunar interior. Thus 
if there were a concentration of Fe or Fe/FeS in the center 
these would be molten, and the Moon would have a molten core. 
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4. THERMAL EVOLUTION OF MARS , VENUS AND MERCURY 

In this section we calculate thermal history models for 
planets Mars, Venus, and Mercury. The Earth has been ex- 
cluded from these calculations. Its evolution history is 
similar to that of Venus and generally is well understood. 

Its present ‘day temperatures are strongly constrained by the 
data, and models can be tailored to meet these requirements. 

New data that became available for Mars, Venus and Mercury 
from recent planetary missions however , justify a re- 
examination of the thermal histories of these planets. 

For the calculations we use the numerical techniques 
described in Section 2 and follow a procedure similar to that 
used for the Moon. Since there are much fewer data for the 
planets, the calculated thermal models are less rigidly con- 
strained than the lunar models. 

4.1 Mars 

The thermal history of Mars is not as strongly constrained 
as for the Moon or Earth. However, new data, primarily from 
the Martian Orbiter Mariner 9, provide information that places 
some conditions on its evolution and internal structure. We 
will now present several theoretically calculated thermal evo- 
lution models from Johnston et al. (1974). 

4.1.1 Structure and Constraints 

The two main pieces of data regarding the composition and 
structure of the interior of Mars are its mass. (M) and moment 
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of inertia factor (I/MR 2 )'. These roughly determine the den- 
sity variation within the planet and may indicate the presence 
(or absence) of a high density core. 

Using values. of and the mean radius as determined from 

Mariner 9 experiments (Lorell, et_ al^. , 19 72) , we calculate 
I/MR 2 = .377 + .001/ indicating a distinct density increase 
with depth. ’This increase could easily be explained by a 
high density core. Recent results from the U.S.S.R. Mars 
probe indicate that Mars maintains a relatively weak dipole 
magnetic field (Dolginov et al . , 1972) giving further credence 
to the idea of a core. 

The presence of a core places a fairly strong constraint 
on the thermal evolution. For core formation to take place, 
temperatures within the planet must have been nearly every- 
where above a core material melting curve at one time. If 
Mars evolved in a manner similar to that of the Earth and Moon, 
the differentiation of core and mantle most likely took place 
relatively shortly after origin. It thus seems to be the case 
with Mars, volcanism and therefore differentiation seeming to 
have been active since early in the Martian history {Carr, 1973). 

We assume that Mars has retained sulfur during its con- 
densation and thus has an Fe-FeS core (Lewis, 1972; Anderson, 
1972) and that potassium is partitioned into the sulphide phase 
(Lewis, 1971; Goettel, 1972) thus providing a heat source in 
the core by the decay of K^. This model leads to early 
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differentiation of Mars and early core formation without 

requiring high initial temperatures necessary for the reduction 
of Fe or Si. This is possible due to the low and relatively 
pressure independent Fe-FeS eutectic at 990 °C (Brett and Bell, 
1969). 

Much new data concerning the evolution of Mars was obtained 
from the Mariner photographs. The general picture was a 
planet much more dynamic than was previously thought. From 
the geologic interpretation of these photographs (McCauley 
et al> , 1972; Carr, 1973; Carr et al . , 1973; Masursky, 19 7.3), 
generally four regions of prominence, stand out: an ancient . 

cratered terrain, younger volcanic plains, large volcanic con- 
structs such as shield volcanoes, domes and craters, and ex- 
tensive sedimentary deposits. 

The Martian evolution inferred from the above features 
can be summarized as follows (Masursky, 1973) : 

1. Roughly 4.6 b.y. ago Mars accreted. An intense bom- 
bardment of the surface took place at, or somewhat later than, 
the final stages resulting in the heavily cratered, oldest 
terrain. 

2. Volcanism, and probably the differentiation of the 
core, followed shortly thereafter., 

3. Intense volcanic and tectonic activity characterized 
by basalt plains, uplifted, faulted and eventually eroded, 
apparently the next stage of Martian evolution. It was in 


was 
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this stage that the volcanic shields such as Nix Olympica 
began to develop. 

4. The next stage featured large scale tectonic activity 
and the completion of the shield volcanoes. 

5. The final stage of evolution, to the present day, 
shows evidence of moderate tectonic and volcanic activity and 
wind and possibly water erosion and deposition. 

Although no absolute ages are available for these stages 
of Martian evolution, they are valuable in constraining thermal 
history models. We know that an initial crust must have dif- 
ferentiated and the core formation must have occurred early in 
the Martian history, with volcanism, implying mantle melting, 
beginning shortly after and proceeding to the present time. 

4.1.2 Heat Sources and Input Parameters 

The thermal evolution models of Johnston et al. , (1974) 
are calculated in the manner described previously with the ex- 
ception that the formation of an Fe-FeS core and the differen- 
tiation of K 4 ^ into the core is specifically allowed. 

As before, the most important heat source in the thermal 
history calculations is the heat generated by the decay of the 
long half-life radioactive isotopes u 238 , u 235 , Th 232 , and K 40 . 
No surface measurements of the abundances of these isotopes 
have been made for Mars so we are forced to look at the only 
available data: from the Earth, Moon, and meteorites. The 
radioactive abundances for the Moon and some meteorites were 
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discussed in Section 2.1.2. The use of lunar, howarditic or 
eucritic compositions in Mars, however, results in present day 
temperatures too hot for any reasonable model. 

The "chondritic coincidence" for the Earth has long been 
recognized (Urey, 1955) giving U * 11 ppb, but this requires 
the depletion of potassium in the crust relative to chondrites. 
Although putting the potassium in an Fe-FeS core can eliminate 
this problem, if all the excess potassium is put in the core, 
this increases the . temperature significantly. In the other 
extreme, assuming no potassium in the core and K/U = 10,000, 
MacDonald (1964) recalculated , the absolute uranium abundance 
for the Earth’s mantle constraining . the heat production under 
steady state conditions to match the average terrestrial heat 
flow. His uranium concentration, reduced to a homogeneous 
planet, is about 31 ppb. One might, therefore, place bounds 
on the Earth ’ s absolute abundances with U from 11 . to 31 ppb, 
and K/U from 80,000 to 10,000. 

Recent models of the chemistry involved in the condensa- 
tion of the solar nebula (Grossman, 1973; Grossman and Larimer, 
1974) indicate that a. number of refractory trace elements, 
including the rare earth’s, will crystallize along with Ca-Al 
rich assemblages at about 1500 °K. This is before the conden- 
sation of magnesium silicates and metallic iron. Potassium 
and other volatiles will condense at lower temperatures (between 
1100°K and 1200°K). Since these temperatures are higher than 
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the condensation temperature for Mars (%450°K) , we might expect 
that the planet retained solar proportions of U, Th, and K. 

This amounts to a bulk uranium content of 15 ppb with K/U = 
50,000 and Th/U =4. We have assumed, however, a bulk radio- 
active abundance equivalent to MacDonald's earth model, i.e. 

U = 31 ppb and K/U = 10,000. It is important to note, though, 
that the total heat production of these abundances is the same 
as that obtained with solar proportions. The only difference 
in the thermal modelling arises in the percentage of potassium 
allowed to enter the suphide phase of the core. For solar 
abundances, this necessarily requires less potassium than is 
assumed in this paper (to be discussed later) but since little 
is known quantitatively about the partitioning of potassium 
and given uncertainties in solar abundances, we believe this to 
be a minor effect. 

The thermal conductivity of the Martian mantle is assumed 
to be temperature dependent and described by the empirical 
equation derived by Schatz and Simmons (1972). The silicate 
melting curve is taken to be an extrapolated peridotite solidus 
(Green and Ringwood, 1969) . Melting, convection of molten 
material and differentiation of the radioactive isotopes is 
modelled as before. 

The formation of the Fe-FeS core begins when temperatures 
in some portion of Mars exceed the Fe-FeS eutectic of 990°C 
and is completed when the temperatures for the bulk of the 



45 . 


planet are above the liquidus for an assumed composition of 
iron and sulfur. In the calculations, core formation is 
modeled as both a continuous and discrete event. When the 

temperatures exceed the eutectic at a grid point, a small 

/ 

amount of excess heat, equivalent to the mass fraction of core 
material, is diverted to the next lower grid point and con- 
verted to its temperature equivalent. Thus, the core material 
effectively melts its way into the interior. Also, a sped-., 
fied percentage of the available potassium is released from 
the silicates in one discrete step and allowed to enter the 
core. Internal heating is also slightly enhanced by the gravi- 
tational settling of the molten metal. The temperature field 

... * 

is reset to an adiabatic gradient when the internal tempera- 
ture rises above the Fe-FeS liquidus where it is believed 
sufficient mass transfer has occurred to warrant this assump- 
tion. The gravitational separation of the core is irreversible 
and once formed, the core does not react with the mantle. 

We have assumed that the Fe-FeS eutectic and liquidus 
are unaffected by high pressure phase changes of FeS, and 
while this is apparently not the case (Usselman, 1974) , we 
maintain that the thermal models and conclusions presented in 
this paper remain essentially unchanged. 

Parameters used in the thermal calculations may be found 
in Table 1. 


i 

t 

\ 
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4.1.3 Thermal Models 

Figure 15 from Johnston et al. , (1974) shows the effect 

of varying uranium concentrations- In the figure, the present 
day temperatures in a uniform, initially cold (T = 0°C every- 
where) Mars, is shown as a function of the present day concen- 
tration of uranium. The terrestrial value of K/U = 10,000 is 
assumed. Also shown are the Fe-FeS eutectic, the Fe-FeS 
liquidus for 85 weight percent Fe and the extrapolated peri- 
dotite solidus but no melting or differentiation is allowed. 

For u 1 37 PP b ' the temperatures are below the solidus. Clearly, 
these models do not satisfy the boundary condition that melting 
and differentiation of the silicates has occurred. Also note 
that for U = 20 ppb, the temperatures do not even exceed the 
Fe-FeS liquidus, prohibiting the completion of a core. To 
satisfy the boundary conditions, with the assumed radioactive 
heat production, higher initial temperatures must be used. 

A thermal model that satisfies the major constraints on 
the evolution of Mars is shown in Figure 15 (model "B” of 
Johnston et aX . ) The initial temperature is calculated from 
the accretion model with an accretion time, t, of 10 5 years 
and a base temperature of 500°C. The core begins to form 
immediately and is completely formed by 1.0 billion years. 

Dry silicates begin to melt at 2 billion years. This zone of 
partial melt expands and works its way into the interior. The 
surface cools slowly by conduction .so that at the present time 
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the dry solidus is reached at a depth of 1800 km. The calcu- 
lated heat flow for this model is: 43 erg/cm 2 -sec , between 
the average measured heat flows for the Earth and Moon. The 
core is presently molten. 

The evolution of Mars and the extent of melting is shown 
schematically for this model in Figure 17. Temperature and 
melting are plotted as functions of depth and time since ori- 
gin. The light shading above the dry silicate melt shows 
regions where the introduction of a small partial pressure of 
water would result in a partial melt. Two episodes of Martian 
evolution are apparent from the figure. The first is the 
early differentiation of a crust due to partial melting in the 
near surface regions and the formation of an Fe-FeS core. 

The formation of the crust is crucial in the understanding of 
the surface geology. Early volcanism, evident from Mariner 
photographs and the age of many regions inferred from crater 
densities, implies the formation and maintenance of a crust 
early in the Martian history. This is followed by upper mantle 
melting and differentiation probably resulting in extensive 
surface volcanism. The region of mantle partial melt pro- 
gresses towards the interior as the surface slowly cools by 
conduction. 

The presence of water and its relation to present day 
tectonic activity is discussed in Johnston et al. (1974). it 
was found that even for models such as the one discussed in 
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this paper, where the anhydrous solidus is not reached until a 
depth of 1800 km, partial melting can occur in the upper mantle. 
It is therefore concluded that one could expect moderate tec- 
tonic activity at the present time. A 200 km thick lithosphere 
obtained for this model is consistent with the depth required 
by hydrostatics to elevate magma to the summit of volcanoes 
like Nix Olympica, about 23 km above the mean Martian sphere, 
and with the gravitational roughness of the planet. 

Although the above temperature model is not a definitive 
statement about the present thermal state of Mars, it is ade- 
quate in the sense that all of the specified boundary con- 
ditions on thermal evolution are met. Hopefully, data from 
the future missions will provide additional constraints and 
thus, thermal models will be better specified. 

4.2 Venus 

Venus is perhaps the most unconstrained planet as far as 
theoretically calculated thermal models are concerned. It can 
be supposed, however, that from size considerations alone, the 
evolution of Venus might be similar to the Earth. There are 
several important differences, though. Condensation models 
of the solar nebula imply that Venus retained little or no 
sulfur, thus eliminating the possibility of FeS or potassium 
in a core (Lewis, 1972). Therefore, proceeding on the bias 
of early differentiation we must require high initial tempera- 
tures in order to form a core within the first billion or so 
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years. If one takes the moment of inertia factor for Venus 
to be equivalent to the Earth's (a rather strong assumption) 
this gives a core of radius about 2900 km for a composition of an 
Fe-Ni alloy. Lewis (1972) also suggests that this core is 
surrounded by a massive mantle of Fe TI -free magnesium silicates. 

A crust, similar in composition to the Earth's, might also be 
expected. 

Radar maps of the Venusian surface indicate relief 
of several kilometers (Campbell et al^. , 1972). Consider- 
ing the density and velocity of the atmosphere we be- 
lieve these features to be subject to extensive erosion 
and are therefore young and probably of tectonic origin. David- 
son and Anderson (1967) have proposed that the rate of volcan- 
ism and tectonic activity is greater on Venus than on the Earth. 
Weertman (1967), however, believes this not to be the case. 
Resolution of this conflict might be found from thermal models. 
Certainly, though,, if the stresses associated with the non- 
hydrostactic state of Venus are not supported statically, we 
would have to expect convective motions in the upper mantle. 

We now present a theoretical thermal model for Venus 
(Figure 18). The computational parameters for this model may 
be found in Table 1. The. initial temperature profile is cal- 
culated from equations (1) and (2) with an accretion time of 
.25 x 105 years- taken to insure core differentiation within 
1.0 b.y. The effect of adiabatic compression on the initial 
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temperature is also included. The total heat production due 
to radioactive decay is essentially chondritic, having assumed 
uranium and potassium abundances equivalent to MacDonald's 
(1964) Earth mantle values, i.e. U = 30 ppb and K/U = 10,000 
at the present day. This is consistent with theories on the 
condensation of the solar nebula (Grossman and Larimer, 1974). 
The bulk uranium concentration for Venus might actually be 
slightly higher than the Earth's, due to the fact it has not 
retained sulfur and therefore has a lower mean atomic weight 
(Lewis, 1972). The thermal conductivity of the undifferen- 
tiated planet and of the differentiated mantle is from the 
Schatz and Simmons (1972) model. 

Core formation is initiated when the temperatures within 
the interior exceed the Fe melting curve calculated for a homo- 
geneous planet following Higgins and Kennedy (1971) . At this 
point (about 1.0 b.y.) it is assumed that the rate of inflow 
of molten iron to the core is sufficient enough so that the 
gravitational energy due to core separation (Birch, 1965; 

Flasar and Birch, 1973) allows for immediate core formation, 
and temperatures in the interior are raised by the heat equiva- 
lent. We reset the temperature profile to an adiabatic gradi- 
ent after core formation due to the extensive mass transfer 
that has presumably taken place. Radioactive isotopes origin- 
ally in the core region are displaced and distributed homo- 
geneously throughout the mantle. 
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.The temperature rise due to core infall is sufficient to 
initiate mantle melting immediately after core formation. The 
mantle melting curve is derived from Lindemann' s equation which 
can be written in a form so that T m = constant • <j > , where T m 
is the melting temperature and $ is the seismic parameter 
(D.H. Chung, personal communication). Using measured values 
of 0 for the Earth, this equation provides an adequate estimate 
for the mantle melting curve for Venus. Convection of molten 
material and differentiation of the radioactive isotopes are 
modelled as described in Section 2.2.2. 

Present day temperatures are shown in the curve labeled 
4.6 b.y. The iron melting curve shown in Figure 18 is calcu- 
lated for a differentiated planet. Notable features of this 
model are the thin (100 km) lithosphere, partially molten 
upper mantle, solid lower mantle, and molten core. Despite 
the molten core, the absence of a magnetic field for Venus 
can be explained by the lack of heat sources in the core, and 
the slow rotation rate of the planet. Calculated heat flux 
is 98 ergs/cm^-sec. 

Figure 19 shows schematically the evolution of Venus as 
a function of time. Isotherms are in degrees C. As noted be- 
fore, core formation takes place at 1.0 b.y. and present day 
conditions indicate a region of high fractional partial melt 
(shaded area) in the upper mantle (higher than expected in the 
Earth) overlain by a lithosphere about 100 km thick. Although 
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the degree of partial melt may be a function of the high sur- 
face temperature, {presumably due to the atmosphere which may 
not have been in existence throughout the hi.story of Venus) 
we have computed thermal models with the surface temperature 
arbitrarily set at 0°C,and this results in a change of partial 
melt only in the upper few grid points. We therefore believe 
that any stresses associated with hydrostatic disequilibrium 
are supported within Venus by strong convective motions in the 
upper mantle. This leads directly to the possibility of 
plate motion, similar to the Earth and likely at a higher rate, 
and we would expect significant tectonic activity at the pres- 
ent time. 

4 . 3 Mercury 

The knowledge of the planet Mercury has increased signifi- 
cantly since the Mariner 10 mission of April, 1974. Prelimin- 
ary results of television reconnaissance of the planet's sur- 
face (Trask, 1974) show it to be remarkably similar to the 
lunar highlands. Regions resembling lunar mares have been 
observed, implying volcanic activity during Mercury's history. 
The presence of a magnetic field (Behannon et al . , 1974) could 
imply a molten iron core. We propose then, that Mercury, as 
the other terrestrial planets, differentiated early, forming 
a crust similar to the Moon, and an iron core. 

Siegfried and Solomon (1974), using methods similar to 
those described earlier in this paper, have calculated several 
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thermal and density models for Mercury. They point out that 
Mercury, however, is different from the other planets in the 
respect that condensation models of the solar nebula imply 
enrichment in refractory elements such as uranium and thorium 
and retention of little or no potassium and other volatiles. 
Furthermore, the thermal evolution is dominated by the high 
conductivity of iron, which makes up approximately 69% of the 
planet by weight. 

Using appropriate initial temperatures, Siegfried and 
Solomon have shown that Mercury can indeed differentiate and 
form a core. Difficulties arise, however, if one assumes no 
heat source in the core and differentiation takes place earlier 
than about 3.0 b.y., in that the core is presently solid, with 
temperatures several hundred degrees lower than the iron melt- 
ing curve calculated following Higgins and Kennedy (1971) . 

This is in apparent conflict with the presence of a magnetic 
field, if such a field is indeed due to an active core dynamo. 
Furthermore, differentiation at times earlier than 1 to 1.5 b.y, 
seems to be required from the nature of Mercury's surface. 

We have therefore calculated a thermal model for Mercury 
(Figure 20) assuming the presence of some sort of heat source 
in the core. We have chosen this to be, as a matter of conven- 
ience, the decay of K 40 . The initial temperature profile is 
calculated from the accretion model with t = 10 3 years and = 
1127 °C (1400°K) , which provides initial iron melting to a depth 
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of about 900 km. Thus, core formation and mantle melting is 
initiated immediately. Other model parameters such as uranium 
and thorium abundances are taken from Seigfried and Solomon 
(1974) and may be found in Table 1. Mantle conductivity is 
taken to be constant and equivalent to that used for the lunar 
thermal model presented in this paper (K = .45 x 10 6 erg/ 
cm-sec °C) , 

As shown in Figure 20, temperature profiles indicate a 
solid lithosphere 200 km thick at 2 b.y. and about 500 km 
thick at the present time. The core is presently partially 
molten and therefore might be able to sustain convective mo- 
tions required to maintain a magnetic field. The lithosphere 
- s probably thick enough to ensure a low level of tectonic 
activity. We would conclude from this model, however, that the 
regions on Mercury's surface that appear like the lunar mare, 
may actually be younger, tectonic activity having been preva- 
lent on the planet for a longer period of time. We have also 
found that a K 40 concentration of .156 ppm at the present time, 
provides a minimum value required to maintain a partially mol- 
ten core. The calculated heat flux is 60 ergs/cm 2 -sec . 

For comparison we show a thermal model of Mercury in 
Figure 21 that contains no heat source in the core. Other 
model parameters are the same as those described above. In 
this case, although the mantle and core differentiate, the 

entire planet has temperatures below the melting curve at the 
present time. 
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Figure 22 shows the evolution of Mercury as a function of 
time for the model in Figure 20. Given the parameters used in 
this model, we conclude that Mercury could have a molten iron 
core if internal heat sources are assumed, and that present 
day temperatures and thickness of the lithosphere limit the 
tectonic activity to very low levels, if any. 
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5. DISCUSSION AND CONCLUSIONS 


The thermal history models calculated for the Moon, Mars, 
Venus, and Mercury are representative of the whole spectrum of 
sizes and compositions of the terrestrial planets. The 
evolution histories are shown together in Figure 23 for a 
comparative analysis. It is clear from this figure that 
although the initial temperatures are similar for these 
planets, the evolution histories and present day temperatures 
are quite different. The size (radius) of the planet, more 
than any other factor, seems to control its evolution and 
the present day thermal state. For the Moon and Mercury, 
the two smallest bodies, the evolution and differentiation 
took place relatively early in their histories. At the 
present both have thick lithospheres and are mostly solid. 
Venus, on the other hand, seems to be still at the peak of 
its evolution. Mars falls in between. 

, A good illustration of the comparison of evolutionary 
history of the planets can be obtained from the value of 
the total thermal energy of the planet as a function of time 
given by : 


E(t) = ^ v P C p T(t)dV + Q(t)^/M 


where E(t) is energy in ergs per gram, T(t) the absolute 
temperature and Q(t) is heat that goes into fusion, at time t. 
M is the mass and V is the volume of the planet. This is 
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illustrated for the planets examined in this paper in 
Figure 24. 

Several conclusions may be drawn from this figure. 

First, the maximum energy corresponds to the peak of the 

evolutionary process. Second, the rate of cooling is dependent 

primarily on the planet's size and somewhat on composition 

(that is, thermal conductivity). 

For example, the Moon, which is the smallest planet, 

has a peak of thermal energy during the first billion years 

after formation. This corresponds roughly to the period of 

large scale magmatic activity as evidenced by the mare basin 

filling. Afterwards, the energy level decreases rapidly as 

the Moon cools. At the present time we would expect the Moon 

to be relatively inactive tectonically. This is illustrated 

by a simple comparison with the Earth, which appears to be 

at the peak of its evolution at the present. The terrestrial 

lithosphere is about 100 km thick while that of the Moon is 

about 700 km. The earth is tectonically active with seismic 

25 

energy release of about E = 10 ergs/year, while the Moon 

s 

is inactive with E . « 10^ ergs/year (Lammlein et al. , 1974). 
While three- fourths of the Earth's surface is covered with 
rocks of 500 m.y. or younger, the youngest crystalline rock 
found on the Moon is 3.16 b.y. 

Mercury, the next smallest planet, also reaches its peak 
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of evolution early in its history. We would expect, upon 
examining Figure 24 (with thermal energy calculated for 
the model shown in Figure 20) that the bulk of the tectonic 
activity occurred before 2.0 b.y. The sharp increase in 
thermal energy is due to core formation. After 2.0 b.y. the 
energy decreases, showing that the planet is presently cooling 
and is probably dormant. If one does not include heat sources 
in the core, as in Figure 21, the slope of E(t) would be 
greater at the present time and the planet would obviously 
be cooler, heat transfer being dominated by the high 
conductivity of iron. 

Mars, on the other hand, intermediate in size among the 
terrestrial planets, reaches its peak of evolution around 
3 b.y,, maintains a relatively constant level of thermal en~ 
ergy for about a billion years and is now beginning to cool. 
From this we can conclude that, while magmatic and tectonic 
activity has existed on Mars since early in its history, most 
of the volcanism and tectonics probably occurred in the last 
2 b.y. We also expect a moderate level of tectonic activity 
at the present time. 

Venus, in many respects similar to the Earth, differen- 
tiated a core early and has maintained a fairly constant level 
of thermal energy throughout its history and can be considered 
to still be at its peak in evolution. Thus, we can expect 
tectonic processes in Venus, similar to those observed on the 


Earth. 
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In interpreting the above results, one must be careful 
to note that the thermal energies calculated for each planet 
were taken from one, non-unique thermal history model. These 
models, however, are reasonable in that they satisfy most of 
the constraints on the evolution of the planets known at the 
present time. Since they represent a bias towards early 
differentiation of all the terrestrial planets, the thermal 
energy turns out to be a useful tool to examine the effect 
of planetary size on the durations of high evolutionary , and 
therefore tectonic, activity. 

In summary, we conclude the following: 

1. The terrestrial planets were heated in some manner 
during or shortly after accretion. 

2. This initial heating allowed for early core forma- 
tion and/or differentiation. Early differentiation is 
evident from rock ages for the Earth and Moon and from photo- 
graphic interpretation of Mars and Mercury surface features. 

3. For planets with large cores such as the Earth, Venus 
and Mercury, the energy released upon core separation is a 
large portion of the total heat budget of the planet. 

4. The Moon is characterized as a differentiated body 
which evolved early in its history. The models that satisfy 
the constraints define a relatively narrow range for the tem- 
peratures in the outer several hundred kilometers, regardless 
of conduction and convective type calculations. Below 1000 km. 
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the control on temperatures is not as great. Temperatures at 
the center of the Moon may be anywhere between 1000 and 1600°C. 
At these temperatures the Moon could have a partially molten/ 
or if Fe and FeS are differentiated, a totally molten core. 

Acceptable evolution models require extensive differen- 
tiation and upwards concentration of radioactive heat sources 

early in lunar history. This differentiation most likely is 
a product of extensive melting, although differentiation due 
to inhomogeneous accretion cannot be ruled out. A model of 
the Moon which is initially hot, and extensively or totally 
molten can best satisfy the constraints. 

5. Thermal models for Mars, assuming an Fe-FeS core 
composition, indicate formation of the core and the differen- 
tiation of an early crust within the first billion years after 
origin in order to satisfy the given constraints on the thermal 
evolution. These models involve large scale differentiation 

of the mantle silicates in the past 2 billion years. This 
would be accompanied by volcanism and outgassing of volatiles 
at the surface. The inclusion of water in the upper mantle 
system leads to partial melting at a depth of about 200 km or 
greater at the present time. 

6. Venus, while not as strongly constrained as the Moon 
or Mars, is characterized as a planet not unlike the Earth in 
many respects. Core formation is allowed to occur by 1 b.y. 
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and present day temperature profiles indicate a partially 
molten upper mantle overlain by a lithosphere approximately 
100 km thick and a molten Fe-Ni core. 

7. Mercury is taken to be similar to the Moon except for 
its large iron core. Thermal models with an assumed heat 
source in the core yield a planet that has a solid mantle to 
about 500 km depth with a partially molten core. Those with 
no heat sources in the core have present day temperatures well 
below the melting curve for iron, heat transfer being dominated 
by the high thermal conductivity of the iron. 

8. The examination of total thermal energy as a function 
of time for each planet provides a method for determining the 
effect of planetary size on the durations of high tectonic 
activity. The Moon, smallest in size, is presently cool and 
tectonically inactive. Mercury is probably inactive at the 
present time, but its peak of activity lasted longer than the 
Moon. Mars, intermediate in size, reached its peak in evolu- 
tion in the past several billion years but is probably moder- 
ately active at the present time. Venus, having been able 

to retain its thermal energy throughout its history, may today 
have tectonic processes similar to the Earth's. 
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TABLE 1 


PARAMETERS 

USED IN 

THERMAL HISTORY 

CALCULATIONS 



Mercury 

Venus 

Moon 

Mars 

Radius (km) 

2439 

6050 

1738 

3389 

Mean Density' 

5.43 

5.25 

3.34 

3.96 

(Lg/cm 3 ) 

C/Ma 2 

- 

- 

0.394 

0.377 

Surface Temp. (°C) 

77 

327 

-20 

-40 

Uranium, ppb 

44 

31 

60 

31* 

K/U ratio 

* 

10,000 

2000 

10,000* 

Th/lJ ratio 

4 

4 

3.6 

4 

Specific Heat 

1.2 

1.2 

1.2 

1.2 

(sil icates) 
(joules/g°C) 

Heat of fusion 

400 

400 

400 

400 

(silicates) 
( joules/g ) 


*See Text 
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FIGURE CAPTIONS 

Fig. 1 Potassium and uranium concentrations in selected 
meteorites/ lunar samples, and terrestrial rocks. 

Lunar rocks are separated into three categories: 
crystalline rocks represented by filled symbols, 
breccias shown by half-filled symbols, and fines 
shown by empt}' symbols. Points for average chondrites, 
howardites, and eucrites are from Mason (1971). 

Sources for lunar values are given in the text. 

Fig. 2 Th/U ratio versus uranium concentration for selected 
lunar rocks and meteorites. Lunar rocks are. divided 
into crystalline rocks, breccias, and fines with the 
same notation as in Fig. 5 . Points for average 
chondrites, howardites, and eucrites are from Mason 
(1971) . Sources for lunar values are given in the 
text. 

Fig. 3 Thermal conductivity of selected materials. Sources 

of measured data (solid lines) are Cremers (1971) for 

3 

Apollo 11 fines, sample density 1.64 g/cm ; Horai 
et a_l. (1970) for Apollo 11 basalt; Murase and McBirney 
(1970) for synthetic Apollo 11 basalt; Schatz and 
Simmons (1972) for sintered, polycrystalline forsterite 
and for single crystal olivine (Fo n , Fa. ,) and enstatite. 

o 6 1 4 

Theoretical (dashed) curves include those of MacDonald 
(1963) for two different values of the opacity e, and 

i 

i 

l 

» 

! 
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Fig . 


Fig . 


Fig . 6 


one proposed by Schatz and Simmons (1972) for 
polycrystalline olivine of approximate composition 

Fo 90 Fa 10 * The ScI ‘ iatz and Simmons curve, labeled 
S above, and K = constant, were employed in this 
paper . 


4 Lunar electrical conductivity profiles of Sonett 
a 1 « (1972) and Dyal et_ al . (1973) and some 

estimates of present-day temperatures in the Moon 
inferred from these measurements. Included are 
the temperature curve of Sonett et: al. (1971) , an 
interpretation of the Sonett et_ £LU (1971) conductivity 
distribution using the temperature-conductivity rela- 
tion of an olivine with little or no Fe 3+ (Duba et al., 

1972) , and the estimates by Dyal et al . (1972) based 

on their conductivity models. 

Schematic diagram of the structure of the Moon as 
inferred from lunar seismic data (Toksdz et al., 

1973) . 

A summary of early lunar evolution and igneous activity 
at the surface based on Rb-Sr , Ar 39 -Ar 40 , and con- 
cordia ages of returned lunar smaples. Samples from 
different missions are designated by "A" for Apollo 
(A—fl is Apollo 11) and "L" for Luna. Sources for 
this summary are presented in the text. 



Present-day temperature profiles for an initially 
cold moon {0°C everywhere) as a function of present 
day uranium concentration. The solidus of anhydrous 
mare basalt (King wood and Essene, 1970) is also 
shown. 


A model for thermal evolution as a function of time 
in an initially molten Moon. Time, in billions of 
years since lunar origin, is indicated by the number 
adjacent to each temperature profile. On this and 
later figures, the Moon is partially or completely 
molten at those depths where the temperature profile 
lies along the solidus. 

Thermal evolution in a Moon accreted in 100 yr . at 
a base temperature of 800°C with an average present 
day uranium concentration of 60 ppb. Symbols are 
explained in the caption to Fig. 8 . 
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Fig. 10 The effect of convection by solid state creep on 

the thermal evolution of the Moon (see text) . The 
initial parameters for this model are the same as 
for Fig. 8 

Fig. 11 Thermal evolution in a Moon derived from inhomogeneous 
accretion (from Tokcdz and Solomon, 197 3) . The 
initial temperature profile assumes a cold ( 0 ° C) 
accretion over a 1000 year time interval. Primary 
zoning of radioactivity' is assumed (see text) , 

Fig. 12 The maximum effect of a thermal "blanket" layer on 
near surface lunar temperatures 100 m.y. after the 
presumed time of the Imbrium event (600 m.y. after 
lunar origin) . Curve B represents a model with 
. initial conditions the same as Fig. 9 , but wrth 
the introduction of a low conductivity, 10 km layer 
at 600 m.y. Curve A is the near surface temperature 
profile for the model in Fig, 9 , with no blanket 

layer. 

Fig. 13 Range of present day temperature profiles for 
theoretical thermal history calculat ions . The 
upper bound is from Toksdz and Solomon (1973) while 
the lower bound is the convection solution of 
Turcotte et al . (1972). 
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Fig. 14 Thermal evolution of the lunar interior as a 

function of time, based on the temperature calcu- 
lations shown in Fig. 9 . Major episodes of 

differentiation and igneous activity are shown on 
top and are based on lunar sample ages. Isotherms 
are in 0°C. Fine and coarse shading denote regions 
of partial and extensive melting, respectively. 
Zones of moonquake activity and high-attenuation 
of S-waves shown on the figure apply to the 
present. 


Fig. 15 Present day temperatures for an initially cold 

(0°C) Mars as a function of uranium concentration. 
Also shown are the Fe-FeS eutectic and liquidus 
and the extrapolated peridotite solidus used in 
this paper (from Johnston et al . , 1974). 

Fig. 16 Thermal evolution of a Mars accreted in 10 5 years 
with T fa = 500°C. U = 31 ppb with K/U = 10,000 
(from Johnston et al., 1974). 

Fig. 17 The evolution of the Martian interior, as a 

function of time, based on the temperature calcu- 
lations shown in figure 16, Isotherms are in «c. 
The light shading above the dry silicate partial 
melt indicates the region where a small partial 
pressure of water would result in partial melt. 



Fig . 


Fig . 

Fig. 

Fig. 

Fig . 

Fig. 

Fig. 


18 A model for the thermal evolution of Venus. The 

5 

accretion time is .25 x 10 years and the effect 
of congressional heating is included in the initial 
temperature. Input parameters for this model are 
described in the text. 

19 The evolution of Venus as a function of time 
showing core formation and mantle melting (shaded 
region). Isotherms are in °C. 

20 A thermal model of Mercury with a heat source 
assumed in the core. Input parameters for this 
model are described in the text. 

21 The thermal history of Mercury with the same 
parameters as for the model in figure 20 except 
that the core is assumed to be free of heat sources. 

22 The evolution of Mercury as a function of time. 
Shading above the core region indicates mantle 
melting. Isotherms are in °C . 

23 A comparison of the evolution of the planets. 

24 Total thermal energy as a function of time for 
the Moon, Mars, Mercury and Venus. 
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